Superconducting Gap Structure of /€-(BEDT-TTF)2Cu(NCS)2 Probed by Thermal 

Conductivity Tensor 
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The thermal conductivity of organic superconductor k-(BEDT-TTF)2Cu(NCS)2 {Tc=lQA K) has 
been studied in a magnetic field rotating within the 2D superconducting planes with high alignment 
precision. At low temperatures (T < 0.5 K), a clear fourfold symmetry in the angular variation, 
which is characteristic of a d-wave superconducting gap with nodes along the directions rotated 45° 
relative to the b and c axes of the crystal, was resolved. The determined nodal structure is inconsis- 
tent with recent theoretical predictions of superconductivity induced by the antiferromagnetic spin 
fluctuation. 
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Since the discovery of superconductivity in organic ma- 
terials about 2 decades ago, the question of the pair- 
ing symmetry among this class of materials is one of 
the most intriguing problems. In particular, the nature 
of the superconductivity in quasi-2D k-(BEDT-TTF)2X 
salts (k-(ET)2X), where the ion X can, for example, be 
Cu(SCN)2, Cu[N(CN)2]Br or I3, has attracted consider- 
able attention. In these layered organics, Shubnikov-de 
Haas oscillation experiments have established the exis- 
tence of a well-defined Fermi surface (FS), demonstrating 
the Fermi liquid character of the low energy excitation. 
The large enhancement of the effective mass revealed by 
the specific heat as well as magnetic susceptibility mea- 
surements suggests the strong electron correlation effect 
in the normal state. Moreover, it was suggested that su- 
perconductivity occurs in proximity to the antiferromag- 
netic ( AF) ordered state in the phase diagram IT] . Since 
some of these unusual properties suggest analogies with 
high- Tc cuprates ]2] , it was pointed out by many authors 
that the AF spin-fluctuation should play an important 
role for the occurrence of superconductivity g,^ . 

Unconventional superconductivity is characterized by 
a superconducting gap with nodes along certain crystal 
directions. Since the superconducting gap structure is 
intimately related to the pairing interaction, its identi- 
fication is crucial for understanding the pairing mech- 
anism. Although the structure of the superconducting 
order parameter of k-(ET)2X salts has been examined 
by several techniques, it is still controversial as we now 
summarize |l|. Results strongly in favor of unconven- 
tional pairing symmetry came from NMR experiments 
of K-(ET)2Cu[N(CN)2]Br, in which the absence of the 
Hebel-Slichter peak and cubic T-dependence of the spin 
lattice relaxation rate 1/Ti were interpreted as an indica- 
tion of d-wave pairing with line nodes |3,H| . The existence 
of the T-linear term in the thermal conductivity at low 
temperatures of k-(ET)2Cu(NCS)2 also supported the 
presence of line nodes B. However, some of the specific 
heat and penetration depth studies on these materials led 



to conflicting results. For example, recent specific heat 
measurements reported a fully gapped superconductivity 
[0 . Since these measurements rely on the T-dependence 
of the physical quantities, it is more desirable to mea- 
sure the in-plane anisotropy of the gap directly in order 
to probe the gap structure. Very recently, such an at- 
tempt was made by STM B and mm-wave transmission 
H experiments. Although both measurements reported 
the strong modulation of the gap structure, they led to 
completely different conclusions on the node directions; 
the former predicts nodes along the directions rotated 
45° relative to the b and c-axes while the latter predicts 
nodes along the b and c directions. In interpreting these 
experiments, one needs to bear in mind that the STM 
spectrum parallel to the 2D plane can be strongly af- 
fected by the atomic state at the edge. Moreover, an 
alternative interpretation was proposed for the mm-wave 
transmission experiments [10|. Thus, the gap structure 
of «;-(ET)2X salts is far from settled and the situation 
strongly confronts us with the need for a powerful direc- 
tional probe of the superconducting gap. 

During the past few years, the understanding of the 
heat transport in the mixed state of superconductors with 
anisotropic gap has largely progressed |ll| . In particular, 
it was demonstrated both experimentally and theoreti- 
cally that the thermal conductivity is a powerful tool for 
probing the anisotropic gap structure [|l2|-|l8[. Thermal 
conductivity has some advantages, compared to other ex- 
periments. First, it is an unique transport quantity which 
does not vanish in the superconducting state, responding 
to the unpaired quasiparticles (QPs). Second, it is a 
probe of the hulk free from the surface effect. Third and 
most importantly, it is indeed a directional probe, sensi- 
tive to the relative orientation among the thermal flow, 
the magnetic field, and nodal directions of the order pa- 
rameter, as we will discuss in detail later. In fact, a clear 
fourfold modulation of the in-plane thermal conductivity 
K with an in-plane magnetic field which reflects the angu- 
lar position of nodes of d^-2 _y2 symmetry was observed in 



YBa2Cu307-5 | |l^,P^ and 2D heavy fermion supercon- 
ductor CeCoIns [|18|, while such a modulation was absent 
in Nb and the B-phase of UPts with an isotropic gap 
in the basal plane p3Hl9| ]. These fact demonstrate that 
the thermal conductivity tensor can be a relevant probe 
of the superconducting gap structure. In this Letter, 
we have measured the thermal conductivity tensor of k,- 
(ET)2Cu(NCS)2 in magnetic field rotating within the 2D 
superconducting planes. The superconducting gap struc- 
ture was successfully determined by the angular variation 
of K. On the basis of these findings, we discuss the nature 
of the superconductivity of k-(ET)2Cu(NCS)2. 

Single crystals k-(ET)2Cu(NCS)2 were grown by con- 
ventional electrochemical method and their approximate 
sizes are 2x1x0. Imm^. The thermal conductivity was 
measured by the steady-stated method with one heater 
and two Ru02 thermometers. The heat current q was 
applied along the 6-direction. The upper inset of Fig. 1 
shows the FS. In the present measurements, it is very im- 
portant to rotate H within the 2D 6c-planes with very 
high accuracy because a slight field-misalignment pro- 
duces 2D pancake vortices which might act as a strong 
scattering center of the thermal current. Special care 
was therefore taken to keep the perpendicular field due 
to the misalignment less than Hd perpendicular to the 
layers, so that the condition for the absence of pancake 
vortices (lock-in state) was always fulfilled while rotating 
H prt. For this purpose, we used a system with two 
superconducting magnets generating H in two mutually 




FIG. 1. T-dependence of the thermal conductivity in zero 
field. The heat current q was applied along the 6-direction. 
Lower inset: The resistive transition at T^. Upper inset: The 
Fermi surface of k-(ET)2Cu(NCS)2 (solid lines). The Fermi 
surface consists of quasi-lD and 2D hole pocket. The dashed 
lines show the first Brillouin zone with kt and kc axes. The 
thin solid lines show the extended Brillouin zone with k^c and 
ky axes in the similar coordinate style of the high- Tc cuprates. 
The node directions determined in our experiment are also 
shown. 



orthogonal directions and a ^He cryostat equipped on a 
mechanical rotating stage with a minimum step of 1/500 
degree at the top of the Dewar. Computer-controlling 
two magnets and the rotating stage, we were able to ro- 
tate H continuously within the 2D planes with a mis- 
alignment of less than 0.01 degree from the plane, which 
we confirmed by the simultaneous measurement of the 
resistivity. 

We first discuss the T- and H- dependence of k. The 
observed T- and H- dependence were very similar to the 
results of ReL Q . Figure 1 depicts the T-dependence of 
K. Upon entering the superconducting state, k exhibits 
a kink and rises to the maximum value just below Tc- 
As discussed in detail in Ref. |g], the enhancement of k 
below Tf, reflects the increase of the phonon mean free 
path by the electron condensation, which is so because 
the phonon thermal conductivity k^^ well dominates over 
the electronic thermal conductivity k*^' near Tc- Figures 2 
(a) and (b) depict the i/-dependence of k in perpendic- 
ular {HI. 6c— plane) and parallel {H\\ 5c— plane) field, 
respectively. In perpendicular field, k(H) shows a mono- 
tonic decrease up to Hc2 above 1.6 K, which can be at- 
tributed to the suppression of the phonon mean free path 
by the introduction of the vortices. Below 1.6 K, k{H) 
exhibits a dip below Hc2- The minimum of k(H) appears 
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FIG. 2. _ff-dependenceof the in-plane thermal conductivity 
(a) in perpendicular and (b) in parallel field {H \\ c) at low 
temperatures. Deviation from the horizontal line shown by 
arrows marks Hc2 ■ Inset: //-dependence of Ak*^ /«:(0) in 
parallel field. For details, see the text. 



as a result of a competition between k^'^ which always 
decreases with H and k*^' which increases steeply near 
Hc2- Then the magnitude of the increase of k{H) be- 
low Hc2 provides a lower limit of the electronic contri- 
bution. As seen in Fig. 2(a), the electronic contribution 
grows rapidly below 0.7 K; k^'/k„ is roughly estimated 
to be >5% at 0.7 K and >15% at 0.42 K, where < 
and Kn are the electronic and total thermal conductivity 
in the normal state above Hc2, respectively. This dra- 
matic increase of K^/Kn is caused by k^'* which decreases 
much faster than k*^' with decreasing T. In parallel field 
with much higher Hc2 (^30 T), k{H) decreases mono- 
tonically at all temperatures. While k{H)/k{0) shows a 
similar if-dependence at 0.71 K and 1.1 K, it deviates 
from this pattern at 0.42 K. Since the electronic con- 
tribution grows rapidly below 0.7 K, this deviation can 
be attributed to k'^K In the inset of Fig. 2(b), we show 

Ak''^{H)/k{0) (= "''"^^i(of' ^°^ ) at 0.4 K, assuming that 

kP^/k{0) has the same iJ-dependcnce. 

We now move on to the angular variation of k as U is 
rotated within the 2D planes. Figures 3 (a)-(c) display 
k{H, 9) as a function of 6 ={q, H) at low temperatures, 
which are measured in rotating 9 after field cooling at 
6* = 0° {H\\ b). The consecutive measurement with an 
inverted rotating direction did not produce any hystere- 
sis in k{H , 9), which demonstrate that the field trapping 
related to the pinning of the Josephson vortices is negli- 
gibly small. At 0.72 K, k{H, 9) shows a minimum at 



9 = 90°. Similar 0-dependence was observed at higher 
temperatures. On the other hand, the angular vari- 
ation changes dramatically at lower temperatures, ex- 
hibiting a double minimum as shown in Figs. 3(b) and 
(c). In all data, as shown by the solid lines in Figs. 3 
(a)-(c), k{H, 9) can be decomposed into three terms 
with different symmetries; k{9) = kq + '*2e + ^48 where 
kq is a ^-independent term, and K2e — C2e cos 29 and 
K4e — C'iB cos A9 are terms with two and fourfold sym- 
metry with respect to the in-plane rotation, respectively. 
The term K20, which has a minimum at HJ-q, appears 
as a result of the difference of the scattering rate for 
QPs and phonons traveling parallel to the vortex and 
for those moving in the perpendicular direction. Since 
a large twofold symmetry is observed even above 0.7 K 
where k^^^ dominates, K20 is mainly phononic in origin. 
In what follows, we will address the fourfold symmetry 
which is directly related to the electronic structure. 

Figures 4 (a)-(c) display K461 normalized by k„ after the 
subtraction of the kq- and K2e-terms from the total k. At 
T=0.72 K, the fourfold component is extremely small; 
1^40 I/kk < 0.1%. On the other hand, a clear fourfold 
component with \C4g\/Kn ^ 0.2% is resolved at 0.52 and 
0.43 K. Since the contribution of k"^' grows rapidly below 
0.7 K and occupies a substantial portion of the total k at 
0.4 K, it is natural to consider that the fourfold oscillation 
is purely electronic in origin. Although {C^ol at 0.42 K is 
as small as 0.2% in Kn, it occupies approximately 1.5-2% 
in K^' and occupies a few % in k'^'(O) assuming k^ /ku ~ 
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FIG. 3. (a)-(c)Angular variation of K{H,e) in |/xoH"|=2 T 
for different temperatures. 9 —{q,H). The solid 

lines represent the results of the fitting by the function 
k{H,9) = Co + C20COS26 + Cab cos AO, where Co, C20 and 
Cis are constants. 
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FIG. 4. (a)-(c) The fourfold symmetry K4,e obtained from 
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tails, see the text. 



0.15. We note that this value is one order of magni- 
tude larger than that in Sr2Ru04 with an isotropic gap 
in the plane ||l^. We now address the origin for the 
fourfold symmetry. The most important issue here is 
" Is the observed fourfold symmetry in k^^ a consequence 
of the line nodes perpendicular to the layer?" . We will 
show that the band structure inherent to the crystal is 
very unlikely to be an origin of the fourfold symmetry. 
First of all, it can be shown by the group theoretical 
consideration that the anisotropic term in k due to a 
fourfold distortion of the FS is of second order relative 
to the leading terms, since the thermal conductivity Kxx 
is a second rank tensor |21|]. In addition, the crystal 
structure of k-(ET)2Cu(NCS)2 is monoclinic and FS is 
nearly elliptic; the fourfold distortion of the FS should 
be very small if it exists. Second, the in-plane magne- 
toconductivity Aa{H) = <j{H) — a{H ~ 0) above Tc is 
undetectably small even at 5 T due to the very strong 
two dimensionality. In fact, the upper limit of Acr/cr(0) 
roughly estimated from the warp of the FS perpendicu- 
lar to the plane is less than 10^^ at 2 T. Thus, as far 
as the Wiedemann- Franz law holds, the fourfold oscilla- 
tion of Kei arising from the magnetoconductance should 
be undetectably small. These considerations lead us to 
conclude that the observed fourfold symmetry originates 
from the superconducting gap nodes [E2|. 

In the thermal transport in the superconductors with 
nodes, the dominant effect in a magnetic field is the 
Doppler shift of the delocalized QP energy spectrum, 
which occurs due to the presence of a superfluid flow 
around each vortex, and generates a nonzero QP density 
of states (DOS) at the Fermi level (Volovik effect) [|3). 
While the Doppler shift increases k*^' with H through 
the enhancement of the DOS, it can also lead to a de- 
crease of K*^' through the suppression of impurity scat- 
tering time and Andreev scattering time off the vortices. 
At high temperatures, the latter effect is predominant, 
but eventually gives way to the former at low temper- 
atures, as demonstrated in high- Tc cuprates |Q. Since 
K*^' increases with H as shown in the inset of Fig. 2(b), 
the enhancement of the DOS is the main origin for the 
i/-dependence of k*^' at 0.42 K |25|. In this case, rotating 
H within the 2D-plane gives rise to the fourfold oscilla- 
tion in hf'- associated with the DOS oscillation ||l5| , [l6| . 
This effect arises because the DOS depends sensitively 
on the angle between H and the direction of the nodes 
of the order parameter, because the QPs contribute to 
the DOS when their Doppler-shifted energies exceed the 
local energy gap. The DOS oscillation with fourfold sym- 
metry gives rise to k'^^ which attains its maximum value 
when H is directed to the antinodal directions and turns 
minimum when H is directed along the nodal directions 
(see Fig. 2 in Ref. |16|). According to Rcf. ||l^, the am- 
plitude of the fourfold symmetry in the d-wave supercon- 
ductors arising from the DOS oscillation is roughly es- 



timated as \Cie\/n''KO) = 0.082 



eiia\ - n nao ^g^^^g ln(^32A/7r;ir). 
Here A is the superconducting gap, F is the QP relax- 
ation rate, vp and v'p are the in-plane and out-of-plane 
Fermi velocity, respectively. Using F '--^ 2 x lO^^s"^, 
2A/fcBTc=3.54, UF -^ SxlO^'m/s, andu^ ~ 2.5xl03m/s, 
gives \Cig\/K''\<S) ~3%. Thus the DOS oscillation by 
Volovik effect yields |C4e|/K'^'(0) which is in the same 
order to the data. 

The fourfold symmetry enables us to specify the node 
directions, which is crucial for understanding the pairing 
interaction. «;4e exhibits a maximum when H is applied 
parallel to the b and c axes of the crystal, showing the 
gap nodes along the directions rotated 45° relative to the 
b and c-axes; the nodes are situated near the band gap 
between the ID and 2D bands (see the upper inset of 
Fig. 1). This result is consistent with the STM experi- 
ments Q. We emphasize here that the determined nodal 
structure is inconsistent with the recent theories based on 
the AF fluctuation. If one assumes an AF fluctuation 
scenario, it is natural to expect the nodes to be along 
the b and c directions. This is because the AF ordering 
vectors become parallel to the b-axis, which would pro- 
vide a partial nesting. If we take the same conventions 
for the magnetic Brillouin zone as the high- Tc cuprates 
with dx2_y2 symmetry (see Fig. 1 (c) in Rcf. S), the 
superconducting gap symmetry of «;-(ET)2Cu(NCS)2 is 
dxy It is interesting to note that superconductivity with 
dxy symmetry has been theoretically suggested based on 
the charge fluctuation scenario [^. Our present results 
may bear implications on this issue. 

We finally comment on the recent heat capacity mea- 
surements which report a fully gaped superconductivity 
||7|. In our view, their temperature range (T > Tc/5) 
is not low enough to conclude the exponential behavior 
of the heat capacity; the measurements at temperatures 
less than Tc/ 10 would be required. 

In summary, the thermal conductivity tensor of k- 
(BEDT-TTF)2Cu(NCS)2 was studied in a magnetic field 
rotating within the 2D superconducting planes. The 
observed fourfold oscillation provides a strong evidence 
of d-wave symmetry. From its sign, the node direc- 
tions are successfully specified. These results place 
strong constraints on models that attempt to explain 
the mechanism of the superconductivity of k-(BEDT- 
TTF)2Cu(NCS)2. 
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